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Abstract 

A method is suggested for interpolating between small-variable and large-variable 
asymptotic expansions. The method is based on self-similar approximation theory 
resulting in self-similar root approximants. The latter are more general than the two- 
sided Fade approximants and modified Fade approximants, including these as par¬ 
ticular cases. Being more general, the self-similar root approximants guarantee the 
accuracy that is not worse, and often better, than that of the Fade approximants. 
The advantage of the root approximants is in their unambiguous definition and in 
the possibility of their construction, even when Fade approximants cannot be defined. 
Conditions for the unique definition of the root approximants are formulated. Several 
examples from high-energy physics illustrate the method. 
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1 Introduction 


A very often met problem in high-energy physics is the necessity of constructing an analytical 
expression uniformly describing a function f{x) in the whole interval of its domain, say, in 
[ 0 , oo], when only asymptotic expansions are known for the small-variable limit a; —)■ 0 and 
the large-variable limit x —)■ oo. The variable x can represent, e.g., a coupling constant. 

The standard way of treating such problems is the use of the so-called two-sided, or 
two-point, or multipoint Fade approximants [1]. In many cases, these approximants provide 
a reasonable interpolation between the small-variable and large-variable limits. However, 
this method has some weak points, as discussed in Refs. HHa. One of them is the occur¬ 
rence of defects, such as spurious zeroes and poles. The other problem is the ambiguity in 
choosing one of the approximant variants from the Fade table. The important limitation is 
the requirement for the compatibility of expansions in the small-variable and large-variable 
limits. The latter means the following. The standard situation in many problems is when, in 
the small-variable limit, one has an expansion in integer powers, x'^, while the large-variable 
expansion exhibits the behavior x^, with a noninteger power jS. Since the large-variable 
behavior of a Fade approximant Pm/n is x^~^, this implies that the integer power M — N 
is not compatible with the noninteger (3. To overcome the problem of incompatibility. Baker 
and Gammel |1] suggested to use the fractional powers of Fade approximants choosing 

the power 7 so that (M — N)'y = f3. The simplest case of the Baker-Gammel method is the 
polynomial approximant P^jiq in a fractional power 7 = / 3 / M . The Baker-Gammel method 
allows one to correctly represent the leading term of the large-variable behavior, although 
the subleading terms not always can be uniquely dehned M. 

In the present paper, we suggest an original method of interpolation between small- 
variable and large-variable expansions. This method allows one to construct an analytical 
expression uniformly approximating the sought function in the whole domain [ 0 , cxo] and 
reproducing both the small-variable and large-variable expansions. The uniqueness condi¬ 
tions are formulated allowing for a unique dehnition of all parameters. The method is more 
general than that of standard Fade approximants as well as the Baker-Gammel method of 
fractional Fade approximants. Therefore, by construction, the accuracy of our method is 
not worse, and often better, than that of the latter methods, with the advantage of being 
uniquely dehned. The approach is illustrated by several examples from high-energy physics. 

2 Method of self-similar interpolation 

Let us be interested in a physical quantity represented by a real function f[x) of a real 
variable x G [0, cxd]. However, the explicit form of this function is not known, since it is 
dehned by complicated equations allowing only for deriving asymptotic expansions in the 
vicinity of two ends, where x —)■ 0 and x ^ 00 . 

For instance, in the small-variable limit, we have 

/(x) ~ /fc(x) (a: ^ 0 ) , ( 1 ) 

with the series 

fk{x) = /o(x) -h ^Onx'^j , ( 2 ) 
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where foix) is a known function. In many cases, the latter enjoys the form 


fo{x) = Ax^ (- 4 ^ 0 ), 


with a being any real number. 

And in the large-variable limit, we can get 

f{x)^f^^\x) (x^oo), 


with the series 


f^P\x) = J2bnX^' 

n=l 


where the powers /?„ are real numbers arranged in the descending order 


l3n+l < (3n (n = 1,2,... ,p - 1) . 


( 3 ) 

( 4 ) 

( 5 ) 


The values of the powers /3„ can be any, either integer or fractional. 

In what follows, it is convenient to deal with the reduced function f{x)/fo{x), normalized 
so that in the small-variable limit. 


fix) _ fkjx) 

fo{x) fo{x) 

we would have the simple asymptotic form 


(x ^ 0) , 


fkjx) 

foix) 


k 

1 - I - . 

n=l 


The interpolation problem consists in constructing such a representation for the sought 
function f{x) that would reproduce the small-variable, as well as large-variable expansions 
(2) and (5), providing an accurate approximation for the whole domain [0, oo]. 

It is worth stressing that this formulation of interpolation problem is rather general. In 
the majority of cases, practically all realistic problems can be reduced to this representation 
employing a change of variables. Also, the small-variable and large-variable limits are con¬ 
ditional, since it is always possible to interchange them by introducing the variable t = 1/x, 
or more generally, t = 1/x^, with a positive /i. 

The method, we suggest here, is based on self-similar approximation theory IMI] com¬ 
bining the ideas of optimized perturbation theory, optimal control theory, dynamical theory, 
and renormalization-group approach. The main ideas of the theory are as follows. The trans¬ 
fer from an approximate form fk to the form fk+i is represented as a motion of a dynamical 
system in discrete time, whose role is played by the approximation order fc = 0,1, 2 .... The 
evolution equation of the dynamical system represents a kind of self-similar relation, where 
the name self-similar approximation theory comes from. The trajectory of the dynamical 
system, by construction, is bijective to the sequence {fk}- The dynamical system in discrete 
time, that is, a cascade, can be embedded into a dynamical system in continuous time, that 
is, into a flow. The convergence of a sequence {fk} to its effective limit is equivalent to the 
convergence of the flow trajectory to a hxed point, which, in this way, corresponds to the 
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sought function /. The motion is governed by control functions guaranteeing fast conver¬ 
gence to the hxed point. The stability of the method is characterized by the map multipliers 
of the related dynamical system. We shall not repeat here all this machinery that has been 
expounded in all mathematical details in Refs. and thoroughly described in review 

articles [l2l|T3], but we shall use the results of this approach. 

Using the self-similar approximation theory for the purpose of interpolation between two 
asymptotic expansions, we come HME] to the self-similar root approximant 

^ =((... (1 + A,xr + + ... + AkX^^^ . ( 6 ) 

Jo{x> 

First of all, we see that setting here all powers Uj = ±1, we can obtain different Fade 
approximants. And if the powers rij = ±1, except the leading that is found from the 
leading term of the large-variable expansion, then we get the modihed Fade approximants 
of Baker and Gammel [1]. However, such a choice of the powers rij is too restrictive and 
arbitrary. The parameters of approximants, to be uniquely dehned, have to be prescribed 
by the available small-variable and large-variable expansions. 

In our previous publications, the root approximants (6) were used so that the powers 
Hi and parameters Aj were dehned through the one-sided expansion, say, the large-variable 
expansion, while the other expansion, e.g., small-variable expansion, was not reproduced. 
The attempts to hnd all values of n* and Aj from the small-variable expansion resulted in 
the equations with multiple solutions. Such a nonunique dehnition of the parameters is, of 
course, unsatisfactory. Now, we aim at generalizing the use of the root approximants (6) in 
such a way that would allow us to uniquely dehne all powers Uj and parameters A* and that 
both the small-variable as well as the large-variable expansions be reproduced. 

We may notice that the root approximant (6) can be identically rewritten as 


IM. 

fo{x) 





Bk-i 



Bk2 

^ink—2 



(7) 


with the parameters 


and powers 


A^f 
B, = ^ 


Aj+i 


(j = l,2 ,...,fc-l) 


mj =j + 1- jrij (j = 1, 2,..., A; - 1) . 

On the other hand, the large-variable expansion (5) can be represented in the form 


( 8 ) 

( 9 ) 


f^A(^x) = bix^^ 



P3-P2 


Expanding the above expressions (7) and (10) in powers of l/x and equating the similar 
terms, we hnd that these expansions are uniquely dehned provided that Eq. (3) is valid, 

= J ( 11 ) 
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(a 7^ A) . 


( 12 ) 


the largest power Uk is given by the relation 

kuk = /9i — a 

and the other powers satisfy the equations 

rrij = f3k-j - f3k-j+i {j = 1,2,k - 1) . (13) 

In this way, all powers Uj of the root approximant (6) can be uniquely dehned through the 
uniqueness conditions 

juj = j + 1 - Pk-j + (dk-j+i {j = l,2,...,k-l) . (14) 

It may happen that the number of terms in the small-variable and large-variable asymp¬ 
totic expansions are not the same, k ^ p. Sometimes, just a single term of the large-variable 
expansion is known {p = 1), while several terms of the small-variable expansion are available. 
How then the uniqueness conditions (14) will be changed? 

Let us assume that just the leading term of the large-variable behavior is known: 

f{x) ~ Bx^ {x —)■ cxd) . (15) 

It is easy to notice that the general expansion (5) is reducible to the asymptotic form (15) 
by setting /?„ = /3. Then the uniqueness condition (14) reduces to the equality 

Uj = ^ {k > p = 1) , (16) 

where j = 1,2,..., k — 1, while the leading power reads as 

(*^= 1 . 2 ,...). ( 17 ) 

All parameters Aj are uniquely defined from the small-variable expansion. 

In the general case, we have k terms of the small-variable expansion and p terms of 
the large-variable expansion. Therefore, to satisfy these expansions, the root approximant 
must be of order k + p, possessing k + p parameters Aj, among which k parameters Aj are 
dehned by the accuracy-through-order procedure from the small-variable expansion and the 
remaining p parameters are dehned from the large-variable expansion. But then we also 
need to have k + p equations for determining k + p powers Uj, while only p terms of the 
large-variable expansion are given. How all powers Uj could be found in such a case? 

Fortunately, large-variable expansions practically always enjoy the following nice prop¬ 
erty. The diherence 

Al3j = I3j - ij = l,2,...,k-l) (18) 

between the nearest-neighbor powers is invariant: 

Af3j = Af3 = const {j = 1,2,..., k — 1) . (19) 

In that case, with the leading power being always given by an equation of type (12) and 
with all remaining powers Uj being dehned by the uniqueness condition (14), we now have 

nk+p = ^ , juj = j + 1- A/S {j = l,2,...,k + p-l) . (20) 

K p 

Then both the small-variable as well as the large-variable expansions can be satished, 
uniquely dehning all parameters Aj, with j = 1,2,..., k + p. 

Below, we illustrate the method by several examples from high-energy physics. 
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3 Supersymmetric Yang-Mills circular Wilson loop 


The A/" = 4 supersymmetric Yang-Mills theory, in the limit of large nnmber of colors N and 
strong t’Hooft conpling A = g^N is taken sometimes as a model for hot QCD im. Since 
there exists an exactly calcnlable expression for the SU{N) circular Wilson loop [T81IT9] . it 
is useful to start with this case, for which the accuracy of approximations can be explicitly 
estimated. 

The exact circular Wilson loop is given by 

H^(A) = Y/,(v/A), (21) 

where Ji is a modified Bessel fnnction of the first kind. In the weak-conpling limit, one has 


W(A) ~ ( 1 - A^/^ + z ^ - ^ 


24 


64 


(A —)■ 0) , 


and in the strong coupling limit. 


3 A-/A (A^oo). 


4\/^ 

Introdncing the change of the variables as 

f{x) = W{\{x)) , \ = x^ 

we obtain the weak-conpling limit 

f{x) ~ (^1 - X + ^ + 

and the strong-conpling limit 


(x —)■ 0) 


/(x) ~ h 
According to rule (20), we hnd 


X 


- 3/2 


4\/^ 


X 


- 5/2 


X —)■ CX) . 


m = ^2 = Ug = 77,4 = ng = 1 , ng = - - . 

Then the corresponding root approximant becomes 

/g (x) = e* (l -h Aix -F A2X^ -f- Aax^ -f A4X^ -f- Agx® - 1 - Agx®) 
acqniring the form of a modihed Bade approximant, with 


Ai — 4 , A 2 — 
653 


A 4 = 


96 


A, = 


15 

Y ’ 


As — 


26 

y ’ 


Afi — 


TT 


( 22 ) 

( 23 ) 

( 24 ) 

( 25 ) 

( 26 ) 

( 27 ) 

( 28 ) 
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The maximal deviation of the root approximant /g from the function, corresponding to the 
exact Wilson loop W, is 0.003. This accuracy is the same as that of the Fade approximant 
P 5 / 7 , studied in Ref. [20] and requiring the knowledge of twice more terms of the asymptotic 
expansions. 

As has been mentioned above, small-variable and large-variable expansions can be inter¬ 
changed by a change of variables. Thus, in the present case, we may use the change 


^(t) = mMt)), A = i. 


(29) 


Then, dealing with the variable t and following the general rule, we get the root approximant 


^ 6 (^) — (1 -|- Bit B2t^ -\- B/it^ -\- B^t^ -\- Bgt^^ 


6\-i/4 


(30) 


with the parameters 


Bi = - , 

2 

A - 2 ° 


n 663 

^ 247r2 ’ 

R - 

-05 — —^ , 


B^ = 


104 


^ 37r2 ’ 


TT^ 


Approximants (28) and (30) coincide with each other. 

This example is also interesting demonstrating how the modihed Fade approximants 
naturally arise in our method. That is, it is shown that our method includes the modihed 
Fade approximants as a particular case. Such a reduction, of course, happens not always, but 
rather rarely. The root approximants ( 6 ) enjoy a more general form than Fade approximants, 
because of which they can provide better accuracy. 


4 Planar cusp anomalous dimension in supersymmet¬ 
ric Yang-Mills theory 


In the A/" = 4 supersymmetric Yang-Mills theory, in the limit of large angle, the planar cusp 
anomalous dimension is linear in angle, with a coefficient F( 5 f) that is the cusp anomalous 
dimension of a light-like Wilson loop, which depends only on the coupling g. The weak- 
coupling and strong-coupling expansions [2IH26] are 


F(^)~4c/- —g - 


737r2 


+ 4C"(3) 


9^ {9 ^ 0) 


(31) 


and, respectively. 


^ , 3 In 2 

r(9) - 29 —— (9 -»00) 


(32) 


The corresponding root approximant reads as 


^b{9)=^9^ I + ^ +^ 29 '^'^ +^35'^^ + Afi' 


- 1/10 

^ , (33) 
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where 


Ai = 


2567r2 

189 


Ao — 


133767r^ 

59535 


Aq — 


32 


6751269 


[5409l7r® + 12859560(^(3)] , 


Aa 


256 


15 In 2 

TT 



As = 1024 . 


This form (33) practically coincides with the Fade approximant P 5/6 given in Ref. [20] , 
differing from it only by 1 %. 


5 Spinor mass in heterotic string theory 


The 50(32) spinor mass in heterotic string theory admits [5ll^ a perturbative weak-coupling 
expansion 

M{g) ~ (1 + 0 . 23 ^ 72 ) (g ^ 0) (34) 

and a dual expansion in the limit of strong coupling, 

M{g) ~ ^ 3/4 ^ 0.351(7-1) (77 ^ 00 ) . (35) 


The root approximant for this case is 


with 


M ; (9) = s '-"* ( ( (1 + 4./) A , g ^) + 9139') . 

Til = 1.472 , A 2 = 3.159299 , Tig = 1 . 


(36) 


This form is very close to the Fade approximant P 4/1 calculated in Ref. [5], the maximal 
difference being only 0.5%. 


6 Ground-state energy for Schwinger model 

The Schwinger model [27l[28] is a lattice gauge theory in (1 + 1) dimensions representing 
Euclidean quantum electrodynamics with a Dirac fermion held. It possesses many properties 
in common with QCD, such as conhnement, chiral symmetry breaking, and charge shielding. 
For this reason, it has become a standard test bed for the study of numerical techniques. 

Here we consider the ground state of the model, corresponding to a vector boson of mass 
M(x) as a function of the variable x = m/g, where m is electron mass and g is the coupling 
parameter having the dimension of mass, so that x is dimensionless. The energy is given by 
the relation E = M — 2m. The small-x expansion for the ground-state energy [291132] is 

P(a;)- 0.5642 - 0.219X +0.1907x2 (x ^ 0) . (37) 

In the large-x limit [321135] . we have 

P(x) - 0.6418x-^/^ - - x-^ - 0.25208x-'^/^ {x ^ 00 ) . (38) 

TT 







The corresponding root approximant is 


/ / / 7/6 \ \ 

El{x) = A I + A2 x'^'^ + A^x^j +A4 x'^] + A^x' 


- 1/15 


(39) 


where 


A = 0.5642 , 
As = 4.028571 , 


Ai = 3.109547 , 
A 4 = 1.070477 , 


A 2 = 3.640565 , 
A 5 = 0.144711 . 


The accuracy of this approximant can be compared to data obtained in other calcula¬ 
tions, density matrix renormalization group, Edmrg, by Byrnes et al. [36] and fast moving 
frame estimates, Epmfe by Kroger and Scheu [37|. This energy has also been calculated 
by using variational perturbation theory [38|, which however has been found to be rather 
complicated and having no advantage over other techniques, such as Bade approximants [35] . 
Adam [39] used renorm-ordered perturbation theory, although his results are less accurate 
than Edmrg and Efmfe- In Table 1 , we compare the latter data with our result E^. Also, 
we present the energy Epa found by means of the Bade approximant All results 

are close to each other, being practically the same in the frame of calculational errors. 


7 Prediction of large-variable expansions 

The great advantage of the method of self-similar root approximants, as compared to Bade 
approximants, is that root approximants can predict the correct behavior of sought functions 
at large variables, being based on small-variable expansions. Such a prediction is principally 
impossible by means of Bade approximants, when the small-variable and large-variable ex¬ 
pansions contain incompatible powers m, for example, when the small-variable expansion 
is in integer powers, as in Eq. (2), while the large-variable expansion, as in Eq. (5), is in 
fractional powers. 

To illustrate this basic advantage, let us consider the ground-state energy for the Schwinger 
model, studied in the previous section. Suppose, only the small-variable expansion (37) is 
available and the powers —1/3 and —1 of expansion (38) are known. But no coefficients 
from the large-variable expansion are given. With the two terms of expansion (37), we can 
construct the root approximant 

El{x) = A ((1 + A^xf/^ + A2X^y^''' , (40) 

in which all parameters are found from the small-variable expansion (37): 

A = 0.5642 , Al = 1.746721 , A 2 = 0.458024 . 

By extending this expression to large variables, we get 

E 2 {x) ~ 0.642616a:“^/^ {x —>■ 00 ) . (41) 

As is seen, the value 0.6426 very well approximates the hrst coefficient 0.6418 of the exact 
expansion (38). 
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Moreover, taking into account this prediction (41), we can construct the root approximant 



1/9 


(42) 


where again the parameters A, Ai, and A 2 are dehned by the small-variable expansion (37), 
resulting in 


A = 0.5642 , Til = 2.245784 , A 2 = 1.336080 , TI 3 = 0.309979 . 
At large variables, Eq. (42) yields 


E;(x) ~ 0.6426161-'''“’ - 0.3229851--’ (i -> oo) . 


(43) 


Both coefficients here are close to those of the exact expansion (38). The value of the second 
coefficient 0.322985 has to be compared with l/vr = 0.318310. 

8 Remarks on scheme and scale invariance 

It is worth noting that self-similar approximation theory [MU, used for deriving self-similar 
root approximants, considered in the present paper, employs the ideas of renormalization 
group in a sense that is different from this notion in quantum held theory, although being 
close mathematically. In self-similar approximation theory, the transfer from one approxima¬ 
tion of order k to another, say, /c -T 1, is considered as a motion in discrete time k. Then the 
sequence of approximations can be treated as a cascade. Embedding the cascade into a how 
makes it possible to pass from discrete time k to continuous time t. The how evolution is 
represented by the renormalization-group equation of the type df{x)/dt = v{x), where v{x) 
is a how velocity. Solving the latter equation iteratively leads to self-similar approximants 
for the sought function f{x). 

In quantum held theory, such as QCD or QED, there is the known problem of scheme and 
scale dependence of truncated series. The sought function f{x) is not scheme-scale depen¬ 
dent, while its truncated series, as in Eq. (2), generally, can be dependent on both. Thus, if 
f{x) is an observable quantity that is a function of the coupling parameter x = OLsj'n^ then 
both the coefficients an, as well as the coupling parameter x, can depend on renormaliza¬ 
tion scheme and scale. The dependence of x on scale comes from the renormalization-group 
equation dx/dt = f3{x), where t is a scale shift and f3{x) is a Gell-Mann-Low function. It has 
been shown that Bade approximants reduce the scale dependence, when the Gell-Mann-Low 
function is taken in one-loop approximation [iQlllI], and a generalized approach [421 - I45] has 
been developed for achieving perturbative scale invariance for any number of loops. 

In the cases we have considered, the situation is different from the mentioned quantum- 
held extrapolation problems, since we investigate the interpolation method, where both the 
small and large variable limits of the sought function are given. In our case, there can arise 
the question of invariance with respect to a parameter, if expansion ( 2 ) is not unique, which 
is analogous to the scheme dependence. For example, there can exist different perturbative 
schemes, associated with different values of a parameter t, so that the coefficients a„ = anit) 
in expansion ( 2 ) depend on this parameter, while the sought function /(x) should not depend 
on such parameters. Then the truncated series fk{x,t) also depend on the parameter, as 
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a result of which the root approximant includes the dependence on t. However, 

in the interpolation problem, the small-variable asymptotic form /o(a;) is given, being not 
dependent on auxiliary parameters, as well as the large-variable limit is also assumed to be 
available, say, as the asymptotic behavior (15), where B does not depend on t. This implies 
that the small-variable limit 


lim 


fix) 


lim 

x^O 


fkix,t) 

foix) 


1 


is hxed, together with the large-variable limit 


lim 

X^OO 


fkix,t) 

fix) 


lim 


fkix,t) 

Bx^ 


1 . 


(44) 


(45) 


For instance, if the small-variable behavior is given by Eq. (3) and only the large-variable 
asymptotic form (15) is available, then the root approximant is 


fl{x,t) = Ax'^ . .{1 +Aixf + A 2 X^Y^‘^ + ... + Akx’^^ , (46) 


with Uk dehned in Eq. (17). The dependence on t comes from the quantities An that 
are expressed through the coefficients Qnit) in the process of the accuracy-through-order 
procedure. Then the large-variable limit (45) yields 

Al^^...{Al + Alf" + As) +... + Akj =B. (47) 


Since the right-hand sides in Eqs. (44), (45), or (47) do not depend on t, the root approximant 
is asymptotically t-invariant. 

Moreover, the root approximants, as has been shown, approximate well the sought func¬ 
tion in the whole interval of the variable x G [o, cxd). More precisely, the root approximant 
f^{x,t) uniformly approximates the sought function fix), with the maximal error e, so that 


fkix,t)-fix) 
fix) 


(48) 


It is easy to see that if there are two approximants, related to two different parameters ti 
and t 2 , such that 


fkjxAi) - fix) 

fix) 


< £i , 


fkiX:t 2 ) - fix) 

fix) 


< £2 , 


then the difference between these approximants is described by the inequality 


fkjxAi) - fkix,t2) 

fix) 


< El 62 ■ 


(49) 


In that sense, the root approximants, within the given accuracy, are approximately scheme 
invariant, that is, invariant with respect to the parameter t labelling different perturbative 
schemes. 
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The dependence on a parameter can also occur in the change of the variables x = x{z, t), 
similarly to the scale dependence in held theory. Then it is possible to expand fk{x{z,t)) in 
powers of the new variable 2 :, getting gk{z,t). The corresponding root approximant gl{z,t), 
after the inverse transformation 2 : = z{x,t), resulting in gl{z{x,t),t) is such that, by con¬ 
struction, it uniformly approximates the sought function f{x), so that 

gl{z{x,t),t) - f{x) 

fix) 

within the maximal error e'. This allows us to classify the root approximants as approximately 
scale invariant, within the given accuracy. 

We would like to note it again that the problem of interpolation is different from that 
of extrapolation. For the latter, the problem of scheme and scale invariance is more compli¬ 
cated. We plan to consider this in future publications. 


< e 


(50) 


9 Conclusion 

We have suggested a general and simple method for interpolation between small-variable 
and large-variable asymptotic expansions. The method is based on self-similar approxima¬ 
tion theory, which allows for the construction of approximations, whose form follows from 
extracting the properties of functional similarity between the given expansion orders. Math¬ 
ematical details of this theory can be found in the cited references. The resulting self-similar 
root approximant (6) makes it possible to satisfy the small as well as large variable expansions 
and to uniformly describe the sought function in the whole domain of its dehnition. 

The general form of the root approximant (6) includes as particular cases that of Fade 
approximants, because of which the accuracy of root approximants is not worse than that 
of Fade approximants. But the root approximants enjoy several advantages, as compared to 
Fade approximants. 

First, a root approximant for given orders k and p is uniquely dehned, while Fade approx¬ 
imants for each given orders, allow for multiple representations. Really, for a small-variable 
expansion of order k and large-variable expansion of order p, it is admissible to construct 
the whole table of Fade approximants Pm/n, with different M and N satisfying the equality 
M + N = k + p+l. Also, when one intends to construct diagonal Fade approximants, one 
needs an even number of terms in a small-variable expansion, while root approximants can 
be formed for any number of such terms. 

Second important advantage is that root approximants make it possible to predict large- 
variable expansions, being based on small-variable expansions. Such a prediction by means 
of Fade approximants is principally impossible, when the small-variable and large-variable 
expansions contain incompatible powers. Moreover, in the case of incompatible expansions. 
Fade approximants cannot be dehned at all. 

We have formulated uniqueness conditions allowing us to uniquely dehne all parameters 
of the root approximant (6) from the coefficients of the given asymptotic expansions. The 
root approximants are shown to be approximately scheme and scale invariant. 

The use of the interpolation formula (6) is convenient for the problems of high energy 
physics, when, due to the duality between weak coupling and strong coupling, there exist 
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asymptotic expansions for both these limits. We have illustrated our approach by several 
examples, demonstrating the generality, simplicity, and good accuracy of this method. 
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Table Caption 


Table 1; Ground-state energy of Schwinger model, for the varying dimensionless param¬ 
eter X = m/g, in different approximations: Density matrix renormalization group, Edmrg: 
Fade approximants, EpA, fast moving frame estimates, Efmfe, and self-similar root ap- 
proximant E^. 
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Table 1 


X 

Edmrg 

EpA 

EpMFE 

El 

0.125 

0.540 

0.540 

0.528 

0.540 

0.25 

0.519 

0.520 

0.511 

0.519 

0.5 

0.487 

0.489 

0.489 

0.487 

1 

0.444 

0.447 

0.455 

0.444 

2 

0.398 

0.396 

0.394 

0.392 

4 

0.340 

0.340 

0.339 

0.337 

8 

0.287 

0.286 

0.285 

0.284 

16 

0.238 

0.236 

0.235 

0.235 


Ground-state energy of Schwinger model, for the varying dimensionless parameter x = m/g, 
in different approximations: Density matrix renormalization group, Ejjmrgi Fade approxi- 
mants, EpA, fast moving frame estimates, Epmfe, and self-similar root approximant E^. 
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